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Abstract 



OO ^ 

^D I The procedure for the STAR Barrel Electromagnetic Calorimeter (BEMC) abso- 

lute calibrations, using penetrating charged particle hits (MlP-hits) from physics 
events at RHIC, is presented. Its systematic and statistical errors are evaluated. It 
is shown that, using this technique, the equalization and transfer of the absolute 
^ ' scale from the test beam can be done to a percent level accuracy in a reasonable 

Y ' amount of time for the entire STAR BEMC. MlP-hits would also be an effective 

C^! tool for continuously monitoring the variations of the BEMC tower's gains, virtually 

without interference to STAR's main physics program. The method does not rely 
on simulations for anything other than geometric and some other small corrections, 
and also for estimations of the systematic errors. It directly transfers measured test 
r> ' beam responses to operations at RHIC. 
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Introduction 



The calibration of calorimeters in collider experiments, using external beams, 
is typically either impractical or even impossible. As a result, the calibra- 
tion is usually performed in situ, using several independent complimentary 
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methods to make sure the cahbration constants obtained by different meth- 
ods are mutually consistent within the systematic and statistical uncertainties 
of each methods. In most large collider experiments, it is quite common to 
use radioactive sources, cosmic rays, and calorimeter hits from high statistics 
physics events with clear signatures and hit patterns [1,2]. For electromag- 
netic calorimeters (EMC), examples of the latter include energy-momentum 
matching for electrons as they are measured in the tracking detectors and 
calorimeter, measuring the positions of invariant mass peaks for known res- 
onances and short-lived particles (vr", r], J/ip, Z°, ...), etc. Most of these 
processes can and will eventually be used at various stages of the running 
STAR experiment for precise calibrations of its electromagnetic calorimeters 
in various energy intervals. However, the "direct" calibrations over a wide en- 
ergy range, using the physical processes above with electrons and/or photons 
in final states, will require a rather long time to obtain sufficient statistics 
and to do the appropriate data analysis. Therefore several "indirect" com- 
plimentary methods, which are described in Ref. [3], will be used first to set 
the channels' gains to the few percent level and obtain the initial calibration 
constants immediately after a module installation during the few first days or 
even hours of running RHIC 

In this paper, we will focus on just one of many indirect approaches which is 
currently considered as the base-line method [3] for the STAR Barrel Electro- 
magnetic Calorimiter (BEMC) calibration, equalization and continuous gain 
monitoring. This method relies on measuring the BEMC towers' responses to 
penetrating charged particle hits from physics events at the running RHIC 
and comparing them to the ones obtained in the BEMC towers exposure to 
an external test beam. 



1 STAR Barrel Electromagnetic Calorimeter 



The STAR BEMC, a sampling scintillator-lead calorimeter with the pseudo- 
rapidity coverage of —1 < r] < +1 and 27r in azimuth ip, has been described 
elsewhere [3-5]. It will be installed within the STAR magnet at the radius 
of ~220 cm just outside the Time Projection Chamber (TPC) and Central 
Trigger Barrel (CTB). 

The BEMC consists of 120 modules of the size At] x Aip ~ 1x0.1 each. The 
BEMC module (Fig. 1) is a stack of 21 scintillator megatiles, separated by 
5 mm lead sheets. Each megatile is divided into 40 light-insulated tiles, 20 
in 77 by 2 in (/?, with a separate readout from each tile, using wave-length 
shifting fibers connected to ~3 m long clear fibers. The light from each set of 
21 tiles of different layers, which creates a single projective towers of the size 
ArjxAip ~ 0.05x0.05, is collected at a single phototube (PMT), one per tower. 



Front Plate (Al) 

Presl 





^ 


M 

Separaling Grooves 

---^ D 

u 




_- Fiber groove 


L 
E 


F 


^s 


--WLS Fiber 




\ 




M 
C 




*\0|rtical 




T 






C 
A 


^Clear Fiber 




C 




PMT-eleclronic BOX 


A 






B 

L 
E 


1 


Clear Fitter 

\ 

\ Liaht Miner 

Optica 1 Connector 



Fig. 1. STAR BEMC module view (left), its cross section (center), and light col- 
lection scheme(right). 

totaling in 4,800 channels for the full BEMC coverage. The two front tiles of 
each tower, with a second set of embedded fibers to transport the light to a 
separate photodetector, provide a "Preshower" (PRS) signal, which is used to 
enhance e-7/hadron separation in the BEMC [6]. To compensate ~20% loss 
of light due to a second fiber, the two front megatiles were made 6 mm thick, 
while the thickness of the remaining 19 tiles is 5 mm each. For perpendicular 
tracks, the total BEMC depth is ~18 radiation lengths (~18Xo) and ~0.8 of 
the nuclear absorption length. A position-sensitive two-layer gaseous Shower- 
Maximum Detector (SMD) is placed in front of the 6th megatile at the depth 
of approximately SXq within BEMC stack. 



2 Requirements to the STAR BEMC calibrations 



The most stringent requirements on the equalization of the BEMC towers come 
from the needs of the lowest-level, fast LO-trigger. Well equalized towers are 
those, having equal responses, in terms of digitized signals, to electron-photon 
hits of, for example, equal energy, E^ or equal transverse energy Et = E ■ sin 6, 
where the polar angle 6 = 2 x arctan(e~''). The usual BEMC contribution to 
the fast trigger is to select events with a high local and/or global energy 
deposit^ above the chosen LO-trigger's threshold(s). Along with other factors, 
the "sharpness" of thresholds for this kind of "high-i^-r" triggers depends on 
how well towers' responses have been equalized. On the other hand, due to 
a finite BEMC energy resolution, which also contributes to the widening of 
thresholds, it does not make sense to equalize tower gains much better than 
an intrinsic BEMC energy resolution. In the range of interest for LO's Ejn- 
thresholds in STAR from ~3-5 to 10 GeV [7], the expected intrinsic energy 
resolution is cte/E :i^ 5-10% [3,5]. This means that the equalization of the 



Either E, or Et, or both. 



BEMC towers to the level of ~3-4% is sufficient in the E'r-interval of ~3- 
10 GeV. 

It would also be correct to compare acceptable statistical uncertainties of the 
BEMC towers' relative calibrations to an intrinsic BEMC energy resolution. 
The relative tower calibration, which is the knowledge of its calibration coef- 
ficient relative to some base gain is equivalent, to some extent, to the tower's 
equalization but for the data analysis rather than LO-trigger, because the un- 
certainties in the calibration coefficients will effectively contribute to the final 
BEMC energy resolution. The other difference is that, while the equalization 
of towers' responses is important in a rather limited energy interval of LO- 
trigger thresholds, the requirements to the relative calibrations are relevant 
to the entire i^^-range of interest for STAR from ~0.5-l to ~50-60 GeV. At 
energies below ~10 GeV, the knowledge of the calibration coefficients with 
the same statistical uncertainty of ~3-4%, as for equalization of digitized re- 
sponses, would be sufficient. Simulations with GEANT [8] have shown that, 
at the energies above 20-25 GeV, the intrinsic BEMC energy resolution is 
almost constant and equal to ~3%. With this resolution, for example, the 
width (RMS) of Z°-peak, AMz/Mz, in e"^e~-decay mode, including its natu- 
ral width, will be ~3.5%. For two times worse BEMC energy resolution of 6%, 
the width of Z^-peak would increase by a factor of ~1.5 to AMz/Mz — 5%. 
The summary of this consideration is that the requirement on the accuracy 
of the relative calibration for each single BEMC tower of no worse than 3-4% 
seems to be valid for the entire energy range from ~0.5-l to ~50-60 GeV, 
particularly early in the program. Of course, over time, one can and will do 
better, utilizing directly electron-photon hits. But the point is that relative 
calibration coefficients with the level of uncertainty above will not significantly 
impact the physics program. 

In STAR, the most restrictive requirements for the absolute BEMC calibra- 
tion arise from the measurements of differential cross sections, which fall 
steeply with transverse momentum, P^. Fits to the SPS data [9] for inclu- 
sive direct-7 and 7r° productions at Pt > 10 GeV/c give the dependence 
of da/dP^ oc Pj^ ~ . The ISR measured spectra [10] at lower P^'s from 
~4 GeV/c fall even more sharply, oc Pj^ .To measure these differential 
cross sections with systematic errors of no more than ~ 10-20%, the BEMC 
absolute scale in the region of interest has to be known at the accuracy of 
-1.5-2%. 

The requirements on monitoring the variations of calibration coefficients over 
time are directly related to the considerations above. Tracking variations of a 
mean gain for the entire BEMC or its any large patch{^ has to be done at the 
accuracy of the absolute scale, i.e. at about 1.5-2%. The statistical errors for 
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tracking gain variations for each single tower can be larger, ~3-4% as for the 
relative calibration and equalization. 

It is also appropriate to mention here the requirement for an initial setting of 
PMT gains after a BEMC module installation in STAR. Formally speaking, 
the method of "MlP-calibrations" , which is the main subject of this paper, 
requires only that initial gains be set so that MIP peaksQ from penetrating 
charged particles are sitting somewhere within digitizers' (ADC) ranges, rea- 
sonably far from their lower and upper limits. Then, only in a few hours of 
running RHIC, the towers can be equalized to the accuracy of a few percents. 
However, so as to not fully incapacitate during these few hours the functional- 
ity of those triggers, which rely on the BEMC signals, it would be desirable to 
have the newly-installed module more or less equalized from the very begin- 
ning. The initial equalization at ~10-15% seems both sufficient and practically 
feasible by, for example, measuring the towers' responses to cosmic rays just 
before a new module installation and then transferring these measurements to 
STAR, using the Light-Emitting Diode (LED) technique [3]. 



3 Calibration scheme 



Many charged hadrons (along with small admixtures of electrons and muons) 
will be produced in every collision at RHIC. In the central region covered 
by the STAR BEMC, these are mostly pions. When striking the BEMC, a 
significant fraction, ~30-40% of high energy charged hadrons, do not deposit 
a large amount of energy via nuclear interactions, instead depositing only 
~20-30 MeV of energy in the calorimeter's 21 scintillator layers due largely to 
electromagnetic ionization. In this paper, we will loosely call all these charged 
hadrons "Minimum Ionizing Particles" (MIP), producing "MlP-hits" in the 
BEMC towers, resulting in "MIP peaks" in the signal spectra. For relativistic 
particles, the position of MIP peak is nearly independent of momentum and 
particle species. This, along with the high yield of charged hadrons, makes it 
attractive to explore the feasibility of using high energy MIPs for the equaliza- 
tion, calibration, and continuous monitoring variations of the BEMC towers' 
gains. 

The calibration scheme using MlP-hits has two stages. In the first stage, a 
sample of BEMC modules is exposed to an external beam, for example at the 
Alternating Gradient Synchrotron (AGS) of Brookhaven National Laboratory. 
The composition of the AGS B2-Line negative test beam is mainly 7r~'s with 
some fractions of other hadrons, electrons and muons of a chosen momentum, 
selected in the range from 0.3-0.5 to 7-8 GeV/c. Therefore the ratios of each 
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tower responses to MIP's and electron's hits, M/P/ e-ratios[^, are measured 
simultaneously (or almost simultaneously for electrons and MIPs of different 
momenta). This makes these measurements completely independent of the 
actual PMT's and ADC's gains, possible attenuation and distortion of signals 
in cables, delay lines, etc., i.e. in the equipment which might be necessary for 
the test beam run but won't be present in the STAR detector, and vice versa. 
In the second stage, after the modules have been installed in their places in 
the STAR, and the RHIC accelerator is producing collisions, samples of MIP- 
hits of particle composition and momenta as close as possible to those in the 
test beam are accumulated for each tower, and the positions of the resulting 
MIP peaks are measured. This step essentially completes the procedure of 
transferring beam-test results to STAR. For those towers exposed to the test 
beam, their responses to electron hitsQ can immediately be predicted, using 
the known MJP/ e-ratios that have been measured at the beam-test stage. For 
all other modules, these ratios are expected to be close to those of the tested 
ones, provided key design tolerances are kept at the module manufacturing 
stage. 

It is clear that charged particles of low transverse momenta are useless in 
the MIP calibration process because their deflection in the STAR's magnetic 
field causes them to enter the BEMC at large angles. Only a small fraction of 
these particles, if any, pass through all 21 scintillator tiles of a single tower. 
On the other hand, if a P^-threshold is chosen too high, the useful event 
rate would be too low because of the steep drop of the particle yield as Pt 
increases. Simple estimates suggest MIPs of P^'s just above ~1 GeV/c as a 
good compromise between yield and utility. In the STAR magnetic field of 
0.5 T, the trajectories of no fewer than ~50-60% of charged particles with 
Pr > 1 GeV/c, produced at the primary vertexQ, will pass through all 21 
tiles of a single BEMC tower. An additional important benefit is that the 
momenta of ~l-2 GeV/c are within the range of the AGS test beam. 

The calibration errors of the procedures described above depend on a num- 
ber of factors. In the rest of this paper, we will discuss in details the most 
important of them. 



^ In this definition, "M/P" is the MlP-peak position which is described more pre- 
cisely later in the paper. The notation "e" is used here for the ratio Se{Ee)/Ee, 
where Se is the mean BEMC signal from electrons of energy Eg- Thus, the MlP/e- 
ratio represents the energy of electrons which would generate in the BEMC the 
same mean signal as MIPs. 

® Of the momenta actually used in the test-beam run. 
^ And not interacted strongly in the BEMC. 



Design tolerances on the light yield variations from the scintil- 
lator tiles within a tower 



The calibration procedure above relies on the assumption that the measured 
M/P/ e-ratios from a few modules, exposed to the beam in a limited energy 
interval, could be extended to the entire energy range of interest and also ap- 
plied to the non-beam-tested modules. To extend the calibrations to a wider 
energy range, one has to be sure that towers' responses were linear or, alter- 
natively, the functional dependence^ S'e(-Ee) were known. An application of 
the beam-test calibration data to non-tested modules requires that the tow- 
ers of the same rj in all modules have similar functional dependences Se{Ee) 
to within ~3-4% if the specification on relative towers' calibration is to be 
satisfied. This in turn implies certain requirements on the uniformity of the 
module design and manufacturing}^ 
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Fig. 2. Residual uncertainties of cal- 
ibration coefficients after the tower 
equalization, using MIPs and ^^Co 
7-source, placed at various locations 
in the stack. Solid markers: no tile 
selection. Open markers: only tiles 
with light yields within ±20% from 
the mean value selected at the qual- 
ity assurance stage. 
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The functional dependence of a sampling calorimeter response, in general, and 
its linearity, in particular, quite strongly depend on the uniformity of light 
yields from scintillator layers within a tower. Individual tower variations of 
light yield from tile-to-tile introduce nonlinearities in the energy responses due 
to the development of mean shower depth with energyp^. Therefore, setting 
limits during module production on the light yield variations from layer to 
layer represents the most critical issue for the maintenance of the BEMC 
calibration and equalization within specified limits over a broad energy range. 



^ Not accounting scaling factors. 

^ It is worth mentioning here that unknown nonlinearities of the electronic PMT- 
ADC chains, in the beam-test as well as in STAR, will go directly to the systematic 
errors of MlP-calibrations. However, these nonlinearities are relatively easily mea- 
sured and accounted in the calibrations. 
^^ See, for example, Ref. [11] for details. 



On the basis of simulations with GEANT, which are illustrated in Fig. 2, it 
has been determined that, after towers' equalization with MIPs, the residual 
uncertainties of calibration coefficients won't exceed the design limits if light 
yield variations were kept within gaussian a < 10%. From Fig. 2, it also follows 
that the equahzation, using MIPs, is expected to provide slightly better results, 
compared to another widely used method with ^^Co 7-source [2]. 



Reconstructing MlP-peak positions in the beam tests and in 
STAR at RHIC 



The compositions of MIPs in the beam-test and in STAR at RHIC will cer- 
tainly be different and won't be exactly known. To accumulate sufficient statis- 
tics in a reasonably short time, the momentum range of the selected MIPs in 
STAR cannot be made as narrow as it was in the beam test. Moreover, STAR's 
0.5 T magnetic field will change MJP/ e-ratios compared to those measured at 
the external beam. All these differences require introducing multiple correc- 
tions and constitute a potential source of systematic errors. Due to substantial 
backgrounds under and in the vicinity of MIP peaks, their positions are to be 
determined from some, not so obvious fits to signal spectra. The actual back- 
ground depends on the environment and hit selection. One expects different 
backgrounds in the module beam test and in STAR at RHIC. This requires 
that MlP-peak fits to be, to a certain degree, robust against background vari- 
ations, variations of the total light yields from otherwise identical towers and 
electronics noise. 

To evaluate the corrections and the associated systematic errors, the various 
hits in the BEMC towers have been simulated, using GEANT. Then, the en- 
ergy deposits in the scintillator tiles have been transformed to the PMT signal 
distributions, accounting for the photon and secondary-electron statistics and 
electronics noise. In the simulations, the "light-yield" from each single tile var- 
ied from 2.5 to 3.5 photoelectrons per perpendicular MIP crossing. Electronic 
noise, (Tnoise Varied from to 2.5 ADC channels. PMT gains have been chosen 
so as to have MIP peaksP^ sitting in the interval from ~12th to 18th ADC 
channels above the pedestal. Typical simulated signal distributions are shown 
in Fig. 3, along with an example of experimental histograms from the BEMC 
Beam-Test-'98 [12]. 

It was found that Vavilov distributions [13] quite satisfactorily approximate 
GEANT simulated energy deposits from mono-momentum charged pions with 
hadronic interactions turned off. To account also for the photon and secondary- 
electron statistics in the PMTs and electronics noise, the following function 
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Fig. 3. Fits to the typical PMT signal distributions, using V^ VT'E'A'-based 
functions (left column) and gaussians (right column). Row 1 (top) - simu- 
lated mono-momentum pions without hadronic interactions; Row 2 - simu- 
lated mono-momentum pions with hadronic interactions; Row 3 - experimen- 
tal data from the BEMC Beam-Test-'98 [12]; Row 4 - simulated signals from 
pions with PYTHIA-HIJING's spectrum dN^/dPr oc exp[-PT/(0.3 GeV/c)] at 
Pt > I GeV/c; Row 5 ~ simulated signals from PYTHIA-HIJING's mixture of tt^ , 
K^ and p/p at Pt > 1 GeV/c; Row 6 - same as Row 5, but with 0.5 T STAR 
magnetic field on; Row 7 - same as Row 6, but with the background from neutrals 
equivalent to HIJING's central Au-Au events. 
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was used to fit the respective ADC signal histograms (top left frame of Fig. 3): 



dy exp 






VMip{x) = a- j -l-^^^y(x{y)-K,f3^) 

FED '^^yV) 



where: 



<l>v(A; K, /5^) is Vavilov distribution VAVDEN from CERN 
library [14]; 
y-PED-MIP , - 2^ 

^(y) = — , Mip — + ^'i>('^' /5 ); 

^^(-^ ^') - ^^'^/.'^ %??.'.' - the mean value of a "bare" 



Vavilov distribution $y(A; k, /5 



'^(y) = \J^'pED + o. ' iy ~ PED) is the width of the 
smearing gaussian; 
PED and apEo are respectively the ADC pedestal position and 
width, known from separate measurements; 
MIP is for the MlP-peak position, which is expected 

to be close to the mean value over MlP-peak with 
already subtracted PED. 

Constants k, (3"^ and s have been fixed to the following values:[^ n = 0.3, 
/3^ = 0.98 and s = 0.078. Out of 3 parameters, a, MIP and a, which were to 
be determined from the fits, only one was actually the goal: that was MIP. 

Function <l>v(A; «;, /5^) is quite complex by itself, and when it is convoluted 
with a smearing gaussian, the fitting procedure using Vmip{x) becomes rel- 
atively slow, taking 15-20 seconds per histogram in PAW [15] with currently 
available computers. To make quick MIP estimates, simple gaussian fits with 
3 parameters: MIP, a and (Tmip, have also been evaluated: 



Gmip{x) = a ■ exp 



{x - PED - MlPf 



2 X a^.jjp 



To account for a background from hadronic interactions, etc. under and in 



^^ The other approach, with k, /3^ and s selected individually for each tower's rj- 
position and pion momentum from the best $y-fits to the simulated MIP energy 
deposits, has also been tested. It was found however that, in the limited momentum 
range of interest, this approach does not yield better results compared to the version 
with fixed constants. 
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the vicinity of MIP peaks in the signal distributions, a function with 3 more 
parameters, bo, bi and 62, has been added to Vmip{x) and Gmip{,x) in fits to 
histograms of Rows 2-7 in Fig. 3: 



V{x) = VMipix) + bo ■ exp (61a; + b2x'^) 
G{x) = Gmip{x) + 60 • exp {pix + b2X^) 



(3) 
(4) 



Fits were usually extended from PED + 3 x apEo to PED + 2.5 x MIP. 
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Fig. 4. Comparison of M/P-parameters from the fits to the means over signal 
distributions for pions without hadronic interactions. Upper row - 3-parameter fits 
to histograms from mono-momentum pions without hadronic interactions (Row 1 
of Fig. 3); Lower row - 6-parameter fits to histograms from mono-momentum pions 
with hadronic interactions (Row 2 of Fig. 3). Momentum dependences for means 
are shown with small asterisks and curves. 

In Fig. 4, the results of fits to the signal distributions from mono-momentum 
pions are compared to the respective means over histograms from non-had- 
ronically interacting pions. Both, Vmip{x) and Gmip{x) 3-parameter fits to 
the histograms of Row 1 in Fig. 3 yield Ml P-values close to the means and 
to each other within a fraction of a percent. The MIPs from 6-parameter 
^(a;)-fits to the simulated pion signals with hadronic interaction^^ (Row 2 
in Fig. 3) are also in good agreement with the means, deviating by just about 



^^ These distributions were expected to be close to the ones from the BEMC Beam- 
Test-'98 [12], which are shown in Row 3 of Fig. 3, because the composition of the 
AGS B2 test beam, after electron rejection, using Cherenkov counters, was mainly 
7r~-mesons with small admixture of kaons, antiprotons and muons. 
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±0.5% (RMS). However in most cases, G(x)-fits to the same distributions 
underestimate parameter MIP by approximately 1%. 

The characteristics of MlP-hit distributions in STAR at RHIC depend on the 
hit selection as well as on the characteristics of the selected events. For track 
selection in each event, the target towers are those entered by one and only 
one charged particle. This particle has Pt > Pj^^ — 1-1.5 GeV/c. In the high 
multiplicity environment of, for example, central Au-Au events, the selection 
conditions could be modified by an additional "isolation" requirement of no 
charged hits in the neighboring towers either. After checking that an extrapo- 
lated track of interest potentially crosses all 21 scintillator tiles of a toweif^, 
the tower's signal is added to a respective histogram. All the selections above 
are easily made with the STAR tracking system [4]. 

The simulation, which are shown in Rows 4-7 of Fig. 3, have been done to eval- 
uate corrections and systematic errors of MlP-calibrations in STAR at RHIC. 
The results are presented in Fig. 5. At the first step (Row 4 of Fig. 3), it has 
been assumed that all accumulated MlP-hits were due to pions only, and the 
contributions to the signals of invisible hitsP^ from underlying events have also 
been neglected. In this case, the only factor, affecting the positions and shapes 
of MIP peaks, is the pion's momentum spreadP^. For pion spectra, it would be 
natural to expect MlP-peak positions close to that of mono-momentum pions 
with Pt equal to the mean over the spectrum, < Pt >= P^*+ 0.3 GeV/c. 
These expectations are shown with the dashed lines in Fig. 5. One can ob- 
serve that the results of V"(x)-fits agree within the errors with the expectations, 
while G(x)-fits yield M/P-parameter values lower by the same ~1% as they 
were in the case of mono-momentum pions with hadronic interactions on. 

After including the PYTHIA-HIJING's [16] spectra of kaons and (anti)protons 
in the MIP samples (Row 5 of Fig. 3),[^ M/P-parameters move up by ap- 
proximately 1-2% for P|."* = 1 GeV/c. However, at P™* =1.5 GeV/c, there 
is virtually no MIP shifts beyond the statistical errors. Such a difference 
between P|."* = 1 and 1.5 GeV/c could be due to some complex interplay 
of M/P (P)-dependences for all three species which, due to different masses, 
reach their minima at different P. Turning on the 0.5 T STAR magnetic field 
makes the simulated distributions in Row 6 of Fig. 3 close to the ones, ex- 
pected in real minimum bias proton-proton and low-multiplicity (peripheral) 
nucleus-nucleus collisions at RHIC. The magnetic field moves MlP-peak po- 
sitions up by ~1% more at P^^ = 1 GeV/c, and by somewhat less values at 
P|«* =1.5 GeV/c. 



^^ Keeping in mind, of course, extrapolation uncertainties due to multiple scattering. 
^^ Neutrals and energy leaks from neighbor towers. 
^^ Effect of STAR magnetic field is still ignored here. 
^^ 66% of pions, 16% of kaons, and 18% of p/p. 
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Fig. 5. Results of fits to signal histograms for pions of PYTHIA-HIJING's 
Pr-spectrum (vr^, Spectrum, B = 0); PYTHIA-HIJING's mixture of pions, kaons 
and pp without (pp, B = 0) and with {pp, B = 0.5 T) STAR magnetic field; HIJING 
mixture of MlP-hits with the neutral background from underlying central gold-gold 
events {AuAu, B = 0.5 T). See text for details. 

To evaluate the effects of invisible hits from underlying events, the HIJING- 
GEANT energy deposit spectra from neutral particles in the extreme case of 
central gold-gold events have been simulated for the towers at various rj (Fig. 6, 
left frame). Then, by combining these spectra and MIP signal distributions 
from Row 6 of Fig. 3, the histograms for Au-Au central events (Row 7 of 
Fig. 3) have been generated in the straightforward way, using formula: 



V{x) 



Vuipiy) Vneutral{x - fj) dfj 



(5) 



where V{x) is the probability density for the total signal in a tower, and 
Vuipiy) and Vneutrai{x — y) are probability densities for the contributions 
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from MIPs and from invisible neutral hits, respectively. 
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Fig. 6. Typical histogram for HIJING-GEANT energy deposit in a BEMC tower 
from neutrals in central Au-Au events (left frame) and the respective PMT signal 
distribution in this tower (right frame). The locations of MIP peaks are marked 
with the dashed histograms. See text for details. 



It should be underlined that backgrounds from neutrals and, for example, 
from hadronic interactions of MIPs themselves affect signal spectra differently. 
Hadronic interactions of MIPs add to the spectrum of purely ionization losses 
another, much wider distribution, not shifting in any direction the MIP peak 
itself. Contrarily, the nonzero energy deposits from neutrals are added on top 
of MIP signals, thus shifting the total signal to the right in virtually every 
single event. Actually, large energy deposits from neutrals (larger than MIP 
peak width of ~2-3 MeV) are not so dangerous, because they just reduce the 
useful statistics, kicking the respective events completely out of MIP peaks. 
However, the contributions of small energy deposits from neutrals are not 
distinguishable in the fits (3)-(4) from electronics noise. As a result, they 
make the "effective pedestal" and MIP peak looking wider and, due to unipolar 
nature of this effect, they shift both, effective PED and MIP by about the 
same amount to the right. 

From the consideration above it follows that, if in fits (3)-(4), the character- 
istics of the effective pedestals were used for PED and (Jped rather than of 
the electronic pedestals measured, for example, with no beam in RHIC, then 
these fits should yield better estimates for M/P-parameter even in the high 
multiplicity environment of central Au-Au events. These effective pedestals 
are measurable by accumulating the signal distributions in each tower from 
the same sample of events used for the respective MIP histograms, and with 
the same hit selection criteria except that no charged particles at all entering a 
tower. Then, the effective PED and apEo could be extracted from the simple 
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gaussian fits as shown in the right frame of Fig. 6. Using effective pedestals in 
the fits to histograms in Row 7 of Fig. 3 yielded MJP-parameters for central 
Au-Au events in the remarkably good agreement with pp, B = 0.5 T (see 
Fig. 5)0. 

The other interesting observation from Fig. 5 is that, within <1%, the ratios 
MIP{RHIC)/Mean{P = 1 GeV/c) for the both, pp and Au-Au events are 
virtually independent of either t] or P^"*, and for the practical purpose could 
be fixed at ~1. 045-1. 050[^ (shown in Fig. 5 with the solid horizontal lines). 
Actually, the corrections bring MlP-peak positions in STAR to approximately 
MIP{P ~ 2.5 GeV/c) for mono-momentum pions with no magnetic field. Of 
course, this combined correction coefficient is based on simulations only, and, 
for example, due to different compositions of MlP-hits in real events, it might 
require some additional adjustment. The bottom line, however, is that all 
corrections themselves are small, and even if just crudely evaluated by these 
simulations, the residual uncertainties should be well below the limit accept- 
able for the STAR BEMC calibration errors. 



An important thing to underline is that, in the simulations, all charged par- 
ticles originated from the vertex exactly at the STAR's center. In the real 
experiment, the vertices are spread over a significant collision diamond, and 
as a result, MIP peaks could be shifted somewhat due to such spread. How- 
ever, this effect is almost purely geometrical, and the respective correction can 
easily be made to the signal from each hit before filling the histogram. 




Fig. 7. M/P/e(7^) -ratio measured 
in the BEMC Beam-Test- '98 [12] 
and simulated with GEANT. Solid 
lines are for MlP/e- and dashed 
lines are for M/P/ ey-ratios, where 
ct = Se{Ee)/ET and Et = EeSinO 
(see footnote 5 in Sec. 3 for other 
notations) . 



0.2 O.-l 0.6 0.8 



The MIP{P = 1 GeV/c)/e-ratios measured in the Beam-Test-'98 are shown 



^^ Actually, in most cases of measuring the energy of a single photon or electron, it 

would be correct to calculate the net signal relative to an effective pedestal rather 

than the electronic one. However in other cases of, for example, measuring the mean 

electromagnetic energy in the event, usage of electronic pedestals would be more 

appropriate. 

19 For y(x)-fits; using G(x) suggests the lower values by ~1%. 
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in Fig. 7 [12]. Functionally, the experimental 77- dependence of MlP/e is con- 
sistent with GSTARP^ simulations: 

MIP{P = 1 GeV/c)/e = A x {1 + 0.0567/^)/ sin 6 (6) 



However, the normalization constants A differ from each other by ~12%. 
This difference might partially be due to some residual experimental uncer- 
tainties, but overestimating of MlP/e by the standard GSTAR at ~5-10% 
cannot be excluded either. After applying the combined correction coefficient 
above to the experimental data, the following formula could be written for the 
MIP{RHIC)/e //-dependence in STAR: 

MIP{RHIC)/e = (264 ± 4,^* ± 13.^. MeV) x (1 + O.OSGr^^)/ sin 9 (7) 



For the time being, until the better data from measurements at an external 
beam are available, formula (7) is suggested for use in the STAR BEMC 
calibrations at RHIC 



6 Statistics and calibration time 



From the analysis of Sec. 5 it follows that, using MIPs from physics events 
at RHIC allows an absolute calibration of the BEMC at a percent level of 
systematic uncertainties. To achieve comparable statistical errors, a sufficient 
number of useful MlP-hits must be accumulated in each tower. MlP-peak 
widths in signal distributions are expected at ~20%. Thus, to measure MIP- 
peak positions with a relative statistical error S, the number of useful hits 
should be Nmip — (0.2/(5)^. Due to the deflection in the STAR magnetic 
field, the projected trajectories of only 50-60% of all charged particles at 
Pnp — 1 GeV/c, which enter a tower, will cross all 21 scintillator tiles. Out of 
them, 60-70% will experience hadronic interactions and produce signals out 
of MIP peaks. As a results, the total number of high-P^ entries per tower 
should be Nmi ^ 6.5 x Nmip ^ (0.5/5^. 

At design RHIC luminosity, the limiting factor for the useful event rate will 
be the data acquisition bandwidth rather than physics cross sections. To min- 
imize MlP-hit accumulation time, the STAR Level-3 (L3) tracking [17] will be 
exploited a) to select among all incoming events those with high-Py tracks, 
pointing out to the BEMC towers; b) to reconstruct parameters of these tracks 
and, probably, even c) to fill histograms. The input L3 event rate could be as 



20 



GSTAR is the configured with STAR setup version of GEANT. 
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high as 100 Hz. However, the actual MlP-hit accumulation rate will also de- 
pend on the type of incoming events. Here we will provide the time estimates 
for three event types: a) pp and b) Au-Au minimum-bias events with at least 
one charged particle, detected in the CTB {Nctb > 1) within its acceptance 
of— l<?7<l,0<(/9< 27r, and c) low-multiplicity (peripheral) Au-Au events 
with 1 < NcTB < 100. In these events, according to PYTHIA-HIJING, the 
average numbers of tracks with Pt > 1 GeV/c within the BEMC acceptance 
would be 0.23, 7.25 and 1.05, respectively. Putting these numbers together, we 
get the running times of Table 1, to achieve the indicated statistical accuracies 
on the MlP-peak positions. 



Table 1 



RHIC running time and statistical errors for MlP-calibrations 


Statistical error 


20% 


10% 


5% 


2% 


1% 


Minimum Bias pp 


23 min 


1.5 hr 


6 hr 


38 hr 


6.3 day 


Peripheral Au-Au 


5 min 


20 min 


1.3 hr 


8.3 hr 


1.4 day 


Minimum Bias Au-Au 


1 min 


3 min 


12 min 


1.2 hr 


5 hr 



It is worth noting that the above estimates are for the time required to get 
every BEMC tower calibrated. However, after the towers are equalized to a 
few percents, the monitoring of the mean gain variation for a patch, consisting 
of Ntw towers, will take by a factor of Nt^ less time compared to what was 
necessary to calibrate every single tower with the same statistical uncertainty. 



Conclusion 



To summarize, the study in this paper has shown that, using MlP-hits, the 
equalization and transfer of the absolute scale from the test beam calibrations 
can be done to a percent level accuracy in a reasonable amount of time for 
the entire STAR BEMC. MlP-hits are also an effective tool for continuously 
monitoring the variations of the BEMC tower's gains, virtually without in- 
terference to star's main physics program. This method does not rely on 
simulations for anything other than geometric and some other small correc- 
tions, and also for estimations of the systematic errors: it directly transfers 
measured test beam responses to operations at RHIC. 
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